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We show that in very dense quark matter there must exist metastable domain walls where the axial 
U(l) phase of the color-superconducting condensate changes by 2ir. The decay rate of the domain 
walls is exponentially suppressed and we compute it semiclassically. We give an estimate of the 
critical chemical potential above which our analysis is under theoretical control. 
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Introduction. — Domain walls are common in field the- 
ory Jl]]. They are configurations of fields interpolating 
between two vacua. If these two vacua are distinct the 
domain wall cannot decay. There are, however, theo- 
ries with only a single vacuum, in which, nevertheless, 
domain- wall configurations exist. The most notable ex- 
ample is the theory of N = 1 axion In such theories 
the decay of the domain wall is possible, though, the de- 
cay rate is often suppressed. It is generally believed that 
there are no domain walls of either kind in the Standard 
Model. It was advocated only recently that long-lived 
domain walls may exist in QCD |}| at zero temperature 
and baryon density. The possibility of very unstable walls 
was noticed earlier in Ref. [||. Unfortunately, no theo- 
retical control is possible in this nonperturbative regime. 
Somewhat related but different domain-wall configura- 
tions were also discussed in Ref. || in the context of the 
decay of the metastable vacua possibly created in heavy 
ion collisions. 

In this paper we show that, in the regime of high 
baryon densities, where relevant physics is under theoret- 
ical control, QCD must have domain walls. Across the 
wall, the U(1)a phase of the color-superconducting con- 
densate varies from to 2ir. Thus, the same ground state 
is on both sides of the domain wall and, consequently, the 
domain wall is metastable. Our proof of the existence 
and the long lifetime of such domain walls relies on the 
following three facts: (i) the instanton density is small at 
large chemical potential, suppressing the effect of chiral 
anomaly, (ii) the U(1)a symmetry is spontaneously bro- 
ken, and (iii) the decay constant of the pseudoscalar sin- 
glet boson is large. All these effects are known, from ear- 
lier studies, to occur in the color-superconducting phases 
of QCD. However, as far as we know, their implication 
for the domain walls has not been explored. These do- 
main walls are similar to the walls of Ref. |3| , which may 
exist, but no definite statement can be made in this case. 
Asymptotic freedom of QCD allows us to assert the exis- 
tence of domain walls in the high baryon density regime 
reliably. The properties of the walls can be determined 
by controllable weak-coupling calculations. 

We shall also make a rough estimate of the critical 



chemical potential above which such domain walls must 
appear within weak-coupling instanton calculation. This 
chemical potential is relatively high (though not unrea- 
sonably high). In the model of QCD with two light fla- 
vors, the critical chemical potential is estimated to be 
about 6Aqcd ~ 1 GeV, close to the scale at which in- 
stanton interactions become relevant ||. We also calcu- 
late the semiclassical lifetime of the domain wall and find 
it to be exponentially long. 

Domain walls in two-flavor high- density QCD. — The 
simplest model with high-density domain walls is QCD 
with Nf = 2 massless quark flavors (u and d). This model 
is a rather good approximation to realistic quark matter 
at moderate densities, such as in the neutron star interi- 
ors. We recall that the ground state of this model at high 
baryon densities is the two-flavor color-superconducting 
state |t],|| , characterized by the condensation of diquark 
Cooper pairs. These pairs are antisymmetric in spin 
flavor and color (a, b) indices: 
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The condensates X c and Y c are complex color three- 
vectors. In the ground state, X c and Y c are aligned along 
the same direction in the color space, and they break 
the color SU(3) C group down to SU(2) C . The lengths of 
these vectors are equal \X\ = \Y\ and can be computed 
perturbatively (see below). 

In perturbation theory, there is a degeneracy of the 
ground state with respect to the relative U(l) phase be- 
tween X a and Y a . This is due to the U(1)a symmetry of 
the QCD Lagrangian at the classical level. This fact im- 
plies that the U(1)a symmetry is spontaneously broken 
by the color-superconducting condensate. Since this is a 
global symmetry, its breaking gives rise to a Goldstone 
boson, which we denote by rj since it carries the same 
quantum numbers as the i] boson in vacuum. 

It is possible to construct the field corresponding to rj 
boson explicitly. The following object, 



S = 17* = X a Y a * 



(2) 
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in contrast to X and Y, is a gauge- invariant order pa- 
rameter. Furthermore S carries a nonzero U(1)a charge. 
Indeed, under the U(1)a rotations 



e nsa/2 q , 



the fields (|l|) transform as X 
therefore E - 



-2m 



(3) 

*X, Y -> e ia Y, and 
E. Thus, the color-superconducting 
ground state, in which (E) 7^ 0, breaks the U(1)a sym- 
metry. The Goldstone mode 77 of this symmetry breaking 
is described by the phase tp of E, 

E = |E|e- 4V . (4) 

Under the U(1)a rotation (^), tp transforms as 

tp — > tp + 2a . (5) 

At low energies, the dynamics of the Goldstone mode tp 
is described by an effective Lagrangian, which, to leading 
order in derivatives, must take the following form, 



L = f[{d tpf-u\d i tpf 



(6) 



This Lagrangian @ contains two free parameters: the 
decay constant / of the 77 boson, and its velocity u. In 
general, the velocity u of the r\ boson may be different 
from the speed of light (i.e., unity) since the Lorentz 
invariance is violated by the dense medium. For large 
chemical potentials p 3> Aqcd, the leading perturbative 
values for / and u have been determined by Beane, Be- 
daque and Savage H: 



f 2 = M 
1 8tt 2 ' 



(7) 



In particular, the velocity of the 77 bosons, to this order, 
is equal to the speed of sound. The fact that / ~ p plays 
an important role in our further discussion. 

It is well known that the U(1)a symmetry is not a true 
symmetry of the quantum theory, even when quarks are 
massless. The violation of the U(1)a symmetry is due 
to nonperturbative effects of instantons. Since at large 
chemical potentials the instanton density is suppressed 
(see below), the 77 boson still exists but acquires a finite 
mass. In other words, the anomaly adds a potential en- 
ergy term Vi ns t(</?) to the Lagrangian (^J), 



L = f 2 [(d tp) 2 - u 2 {d^f] - V iast (p) . 



(8) 



The curvature of Vinst around tp — determines the mass 
of the 77. 

A standard symmetry argument determines periodicity 
of Vi ns t(tp). One can formally restore the U(1)a symme- 
try by accompanying (0) by a rotation of the ^-parameter 



Nfa = 6 + 2a. 



(9) 



This symmetry must be preserved in the effective La- 
grangian, so the latter is invariant under (H) and (||). 



This means that the potential Vjnst is a function of the 
variable tp — 9, unchanged under U(1)a- Since we know 
that the physics is periodic in 9 with period 2ir, we can 
conclude that, at the physical value of the theta angle 
9 = 0, Vj ns t is a periodic function of tp with period 2ir. 

Moreover, at large /j, Vi nst can be found from instan- 
ton calculations explicitly. The infrared problem that 
plagues these calculations in vacuum disappears at large 
p; large instantons are suppressed due to Debye screen- 
ing [jlO jll I . As a result, most instantons have small size 
p ~ 0(/j, _1 ) and the dilute instanton gas approximation 
becomes reliable. One-instanton contribution, propor- 
tional to cos((y9 — 9), dominates in Vi ns t- Therefore, 



Vinst (<p) 



-ap 2 A 2 cos tp , 



(10) 



where A is the Bardeen-Cooper-Schrieffer (BCS) gap, 
and a is a dimensionless function of p which will be found 
later. Here we only note that a vanishes in the limit 
p, — * 00. This is an important fact, since it implies that 
the mass of the 77 boson, 



m 



a p 
2 7 



A = 2T:y/aA , 



(11) 



becomes much smaller than the gap A at large p. In 
this case the effective theory (||) is reliable, since meson 
modes other than 77 have energy of order A, i.e., are much 
heavier than 77 and decouple from the dynamics of the 
latter. 

The Lagrangian (||) with the potential (|lTJ) is just the 
sine-Gordon model, in which there exist domain-wall so- 
lutions to the classical equations of motion. The profile 
of the wall parallel to the xy plane is 



tp = 4 arctan e 



mz/u 



(12) 



so the wall interpolates between ip = &t z = — 00 and 
tp = 2n at z = 00. The tension of the domain wall is 

(13) 



a = 8v 2aufpA . 



A good analog of this domain wall is the N = 1 axion 
domain wall, which also interpolates between the same 
vacuum. 

Decay of the domain wall. — It is important to under- 
stand the mechanism of the decay of the wall. It has 
nothing to do with the decay of 77 meson quanta, which 
are due to 77 coupling to photons, ungapped quarks, or 
the gluons of the unbroken SU(2) C subgroup. The domain 
wall is already a local minimum of the energy, and the 
decay of its excitations means only that the fluctuations 
around this minimum, corresponding to deformations of 
the wall, are damped. 

The domain wall is not stable because the same ground 
state is on both of its sides: tp = and tp = 27r are equiva- 
lent. The instability is due to higher energy meson modes 
integrated out and not present in the Lagrangian (H). 
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One can visualize the effect of these modes by consider- 
ing an effective potential which depends on the magni- 
tude |E| as well as on the phase <p of the order parameter 
S. This potential has the shape of a Mexican hat, slightly 
tilted by an angle proportional to a. A similar picture 
is discussed in Ref. [||, except that the tilt of the hat is 
very small in our case. The domain wall is a configura- 
tion that, as a function of the coordinate perpendicular 
to the wall, starts from the global minimum, goes along 
the valley, and returns to the starting point. One can 
continuously deform this configuration into a trivial con- 
stant one by pulling the looplike trajectory over the top 
of the hat. 

As in the case of the axion wall |§,|l2| , this deformation 
has to be done in a finite area of the wall first, thus 
creating a hole. If this hole exceeds the critical size, it 
will expand, destroying the wall. On the rim of the hole 
the magnitude of |E| vanishes. The field configuration 
around the rim is a vortex: on a closed path around 
the rim, <p changes by 2ir. The decay of the wall is a 
quantum tunneling process in which a hole bounded by a 
closed circular vortex line is nucleated. The semiclassical 
probability of this process is 



r 



exp 



16tt v z 



(14) 



where v is the tension of the vortex line in the limit of 
massless boson, m = 0. The factor 1/it in the exponent 
of Eq. (|l4|) is due to the fact that u plays the role of light 
speed for the effective dynamics of the Goldstone boson 
[see Eq. (§]. 

To find r we still need to compute the vortex tension 
v. Since the vortex is a global string, its tension is loga- 
rithmically divergent, 



v = 2iru 2 f In 



I? 



Be 



2ttu 2 / 2 ln(-RA) 



(15) 



where R is a long-distance cutoff to be specified later, 
and i? C orc is the size of the core of the vortex line, which 
is the short-distance cutoff. i? C orc ~ 1/A since A is the 
momentum scale at which the effective Lagrangian de- 
scription breaks down. We are helped by the fact that 
the vortex tension is dominated by the region outside the 
core, so the effective Lagrangian (g) is sufficient for com- 
puting v to the logarithmic accuracy. By using Eqs. Jl3| ) 
and (|l5]), and taking R to be the thickness of the wall, 
we find the decay rate to be 



exp 



f 4 



hi' 5 



1 



3 a fj?A 2 ^/a 



(16) 



Since / ~ /j > A, and a decreases with increasing p,, the 
decay rate is exponentially suppressed at high p. Thus, 
we have shown that (i) domain walls exist in the limit 
of large chemical potentials, and (ii) they are metastable 
with parametrically long lifetime. These conclusions are 



valid in the regime of very large chemical potentials p, 
where our calculations are under control. 

Calculation of the potential. — What happens at 
smaller pi The most interesting possibility is that the 
walls persist down to p = as advocated in Ref. || us- 
ing large--/V c arguments. Another possibility is that, as 
the description based on the Lagrangian (||) breaks down, 
the walls disappear. This happens when the mass of the r/ 
excitation becomes comparable to 2A, the typical energy 
scale for higher mesons [Ol. From Eq. (11), one derives 



the following condition when our effective Lagrangian de- 
scription is under control: 



a(p) < 1/tt 2 



(17) 



We shall now evaluate the function a(p). 

To compute Mnst(^)j we start from the instanton- 
induced effective four-fermion interaction fMUqJl5ll, 



Anst = / dpn (p)[ ^ 2 P 3 



(u R u L ){d R d L 



3 

32 



(u R X a u L )(d R X a d L ) 



■^{uRO-^ u \ a u L )(d R o- piu \ a d L ) 



H.c. 



(18) 



By taking the average of Eq. ( |18[ ) over the superconduct- 
ing state (|l|), one finds Vinst, and confirms that it is pro- 
portional to cos ip as in Eq. (|l0|). In the ground state, 



IXI = IYI = i A{po) 

11 11 2 J (27r)*pg + (|p|-/i)2 + Aa(po) 



, (19) 



where A(p ) is the momentum-dependent BCS gap. Us- 
ing the perturbative result p6| , 



A(po) 



we find 



A cos 



9 



(3- : "" A 



la ^ 



1X1 = 



3 p 2 A 
2V2n g 



A < po < fi , (20) 



(21) 



Averaging Eq. ( Jig ) in the superconducting background, 
we find, after some calculations 



V inst ((p) 



dpn (p)[ ^ttV 



12|A| 2 cos^. (22) 



Using the standard formula for the instanton density at 
finite chemical potential pUEjl 



no(p) 
with 



C 



N 



8tt 2 



2 AT 



p exp 



8tt 2 



(23) 
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Cn 



0.466e 



-1.679JV C 



1.34^ 



(JV c -l)!(JV c -2)! 



we arrive at the final result 



a = 5 x 10 4 In 



A 



QCD 



7 / . N 29/3 

M 



(24) 



(25) 



Thus a — > when /i — > oo, so at sufficiently large p, 
the criterion (|lj]) is satisfied. However, due to the large 
numerical constant in Eq. ( p5| ) the critical p is quite high: 
Merit ~ 6Aqcd ~ 1 GeV. This result should be taken 
with some care due to the uncertainty in the criterion 
(|T7|). However, since a depends quite sensitively on p, it 
is reasonable to expect that our estimate is not far from 
the true value. 

In the estimate above we neglected the contribution 
from large instantons, which arise from the unbroken 
SU(2) C sector of the theory. This sector is governed 
by a pure Yang-Mills theory with the confinement scale 
Aq CD ~ Aexp[— const • p/(gA)] jlTj]. The nonperturba- 
tive contribution of large SU(2) C instantons is of order 
(Aq CD ) 4 . Since Aq CD is exponentially small, this contri- 
bution is negligible compared to that from small SU(3) C 
instantons. 

Inclusion of quark masses does not change the domain 
walls in a substantial way. The mass contribution to the 
potential has been found in Ref. M, 



= —b m u md A 2 cos ip , 



(26) 



where b ~ 1. Eq. ( j2q ) has the same <p dependence as Vinst, 
therefore, in all previous formulas one should replace a 
by a + bm u m d / p 2 . 

Discussion. — It would be interesting to investigate 
possible astrophysical consequences of the high-density 
QCD walls. In particular, one would like to know if such 
walls can be created inside neutron stars. To describe 
the motion of the wall one may need more than just the 
effective Lagrangian (@): the coupling of rj to ungapped 
quarks and SU(2) C gluons could be important. The mov- 
ing wall may radiate quark-hole pairs, gluons, or photons, 
slowing down the collapse of a closed domain wall surface. 

It is possible to generalize our results to the color- 
flavor-locking state of N f = 3 QCD The U(1) A 
symmetry is also spontaneously broken in this case. The 
role of the rj boson is played by the rj' meson, which is 
also light at high densities |l9|| . The instanton-induced 
77' potential has a form similar to ( |Io|) [f2(]|j : 



V inst ((p) = -a ■ (m s /n) p 2 A 2 cos ip , 



(27) 



where the evaluation of dimensionless function a' is very 
much the same as our calculation of a and amounts to 
inserting extra factor m s p into (p2|). Thus one expects 
the domain walls to exist and to be metastable at large 
ji. Due to the mixing between the neutral mesons p9[, 



the 7T° and rj fields are also nontrivial on the wall. The 
domain walls also exist in QCD with large isospin density 
pjj . This case is interesting, since it can be studied by 
a Monte Carlo lattice simulation. 
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